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Abstract. We propose a simple stirring experiment to generate quantized ring 
currents and solitary excitations in Bose-Einstein condensates in a toroidal trap 
geometry. Simulations of the 3D Gross-Pitaevskii equation show that pure ring 
current states can be generated efficiently by adiabatic manipulation of the 
condensate, which can be realized on experimental time scales. This is illustrated 
by simulated generation of a ring current with winding number two. While solitons 
can be generated in quasi-lD tori, we show the even more robust generation of 
hybrid, solitonic vortices (svortices) in a regime of wider confinement. Svortices 
are vortices confined to essentially one-dimensional dynamics, which obey a similar 
phase-offset— velocity relationship as solitons. Marking the transition between 
solitons and vortices, svortices are a distinct class of symmetry-breaking stationary 
and uniformly rotating excited solutions of the 2D and 3D Gross-Pitaevskii 
equation in a toroidal trapping potential. Svortices should be observable in dilute- 
gas experiments. 



PACS numbers: 03.75.Fi, 05.45.Yv, 42.50.-p 



Persistent currents, quantized vortices, and solitons in Bose-Einstein condensates 
(BEC) are high-energy excitations predicted by mean-field theory. Quite recently, both 
solitonic waves |], |2| and singly quantized vortices ||, |J, U have been unambiguously 
identified experimentally in 3D harmonic trap geometries. The next step in increasing 
our understanding of quantized currents and solitons is to study their stability 
properties and measure their intrinsic lifetimes thereby testing the validity of mean- 
field theory in a so-far unexplored regime. The torus presents an ideal geometry 
in this context as solitons can travel without reaching an end of the trap and ring 
currents are stabilized by the ring geometry as opposed to vortices where the core 
can drift Q . Excited states in toroidal traps also have the potential to provide useful 
applications. E.g. Drummond et. al. |Q have suggested constructing a mode-locked 
atom laser using a dark soliton with a nodal plane (also called a black soliton or kink 
configuration) in a ring geometry. They demonstrated the creation of a black soliton 
by stirring a condensate in a ID ring with a paddle at constant angular velocity and 
damping additionally generated excitations by Raman outcoupling. Other simulations 
involving nonlinear excitations in toroidal geometries can be found in references j|, §1 • 

In this letter we propose an adiabatic stirring technique to impose circulation onto 
a BEC in a 3D toroidal geometry. Without any need for external cooling or damping 
mechanisms, adiabatic manipulations allow one to generate pure ring-current states 
with essentially arbitrary circulation. A variation of the stirring technique can produce 
solitons if the transverse trap confinement in the toroidal tube is quasi-lD|. For slightly 
less tight confinement new nonlinear objects are created, which have properties of both 
vortices and solitons. We call these objects solitonic vortices, or svortices, and briefly 
discuss their properties. 

In our simulations, the condensate is stirred with a rotating paddle which can 
easily be realized with a blue-detuned laser-generated light sheet. Being purely optical 
in nature, this method does not rely on specific magnetic field configurations and is 
suitable for single-component as well as spinor condensates. Toroidal trap geometries 
have already been used in early BEC experiments 12 and in the recent JILA vortex 

| We define a trapping geometry quasi-lD when the confinement in two of three spatial dimensions 
is such that the ID Gross-Pitaevskii equation essentially g over ns the dynamics, which is the case for 
transverse confinements of the order of the healing length [hoL . 
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experiments ||, |J. Quasi-ID confinement in a toroidal trap certainly presents a 
challenge for experimentalists but seems realistic after promising results in linear 
geometries |l3| and with the prospect of a tightly confining purely optical toroidal 
trap@. 

We model the dynamics of the BEC with the time-dependent Gross-Pitaevskii or 
nonlinear Schrodinger equation (NLSE), which provides the zero-temperature mean- 
field theory |15| ]. For a BEC of N atoms of mass M, confined in an external potential 
V, the NLSE reads: 

ihdtip = [-(h 2 /2M)V 2 + V + g\i/j\ 2 ]i/j, (1) 

where g = Airft 2 a$N / M is the coupling constant. The s-wave scattering length 
is assumed to be positive, relating to repulsive interparticle interactions. The order 
parameter or condensate wavefunction tp is normalized to J \i[)\ 2 dr = 1 and may be 
written as ^>(r, t) = ^J~p exp i<f>, where the square root of the density Jp~ and the phase 
<j) are real functions of r and t. 

We have solved the time-dependent NLSE (Q) numerically on a grid using 
a pseudo-spectral FFT representation and a 4th-order variable-step Runge-Kutta 
integrator in box boundary conditions. The trap potential V — Vho + V co + VJ> a (i) 
mimics a torus by a 3D harmonic-oscillator potential Vj 10 in the shape of an oblate 
disk pierced by a Gaussian blue-detuned laser beam modelled by the core potential V co 
resembling the trap used in early MIT experiments fll2|] . A time-dependent "paddle" 
potential V pa ,(t) mimicking an ellipsoidal Gaussian light sheet is used to stir the 
condensate in the torus|. 

For a condensate of 10 6 atoms of Na with ao = 2.75nm, our simulations relate to 
a torus of radius rx = 32.4/xm (distance of the potential minimum from the symmetry 
axis of the torus). The healing length £(r) = 1/ ^/8iraoNp(r), which depends on 
the spatially inhomogeneous density p(r), sets the relevant length scales for the size 
of solitons (~ 2£) and the transverse size of vortex filaments (ss 3£) Jl6f . From 
the numerically obtained ground state wavefunction in the torus, we find a value 
of £ w 2.2pm for the size of the healing length at peak density. The transverse 
confinement of the condensate in the toroidal tube is best characterized by the number 
of healing lengths, which we define by the line integral = / c £(r) -1 ds, taken 
across the toroidal tube along a radial line C embedded in the symmetry plane of 
the torus. In other geometries, the most sensible curve of integration C is dictated 
by the confinement. The dimensionless confinement parameter J\fc can be computed 
directly form the numerical wavefunction and does not rely on the Thomas-Fermi 
approximation, which is often used to characterize the extent of a trapped condensate. 
In our simulation the confinement in the radial direction is about the same as in the 
axial direction and amounts to ~ 25. Figure [j] shows the result of the simulation 
of a stirring experiment to generate a ring current with winding number w = 2. The 
ground state of the condensate in the toroidal trap intersected by a narrow ellipsoidal 
Gaussian paddle potential was found by imaginary-time propagation. During the 

§ In detail, the harmonic part of the potential is Vho = M{u> x x 2 + u> y y 2 + oj z z 2 )/2 with uj x = uj y = 
OAlUz ; the core-potential is V co = V c exp{— (x 2 + y 2 ) / (2Ar 2 )}. The paddle potential has the form of 
a Gaussian ellipsoid centred at a radius r p and initially reads V pa (0) = V p exp{— (x — r p ) 2 / (2Ar 2 ) — 
y 2 1 (2Aa 2 )}8(x). During the simulation it is accelerated to a uniform rotation around the z-axis 
and retracted by pulling out to the side. The following parameters lj x = 27r8.63Hz, V c = 18.5nK, 
Ar c = 15.1/im, V p = 6.95nK, Aa p = 3.37/^m, Ar p = 27.0/im, r p = 32.4/im, were used in the 
simulation and are comparable to currently achievable experiments. 
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Figure 1. 3D simulation of ring current generation with w = 2. Parts a) and b) 
show gray-scale coded 2D slices of the density p and the phase <f> respectively in 
the xy plane of the trap. Rows c) and d) show the density and the phase on the 
line of minimal potential parameterized by the angle a. Note that the phase is 
plotted within an arbitrary offset and modulo 2ir. The columns show time slices 
of our simulation. The first column shows the condensate ground state in the 
toroidal trap intersected by a resting stirring paddle. Columns 2 and 3 show how 
a phase gradient builds up while the stirring paddle is accelerated at a constant 
rate between t = and t = 0.27s. Afterwards the condensate is stirred with 
constant angular velocity u) = 27r0.10Hz. At this stage, the density notch goes 
down to about 1% of the maximum density and a phase step of A</> m in has 
built up across the paddle. The fifth column shows how the paddle is retracted 
between t = 0.41s and t = 0.68s. The last column shows the w = 2 ring current 
that has established in the toroidal trap. The transverse confinement in the torus 
is m 25 healing lengths. 



simulation the paddle was accelerated to azimuthal rotation "stirring" the condensate 
and subsequently retracted in the direction of the torus axis. 

In order to elucidate the effect of adiabatic stirring on the condensate 
wavefunction and to interpret the result of the simulation shown in figure [j], we 
start with an idealized 1-D picture, neglecting the transverse dimensions of the torus: 
When blocked by a paddle, the torus can be treated as a ID box with hard-wall 
boundary conditions. The ground state ip„ of the condensate in a ID box is a Jacobi 



elliptic function of type sn 17 , which is real- valued and therefore has a constant 
phase. When the paddle is moving, the wavefunction acquires a phase ramp. This 
can be understood by considering the Galilean invariance of the NLSE: The boosted 
ground state ip g (x — vt,t) exp[iM/h(vx — ^v 2 t)] satisfies the boundary conditions in 
the moving frame and is the "ground state" of the translated box. A phase gradient 
in the wavefunction is equivalent to a supercurrent with velocity v — Kd x <j)/M. A 
state like this can either be reached by an adiabatic transition between the resting 
and the moving box through accelerating the paddle slowly or one can prepare the 
condensate right from the start in a torus with a rotating paddle. Removal of the 
stirring paddle reconnects both ends of the "ID box" to a ring and can be modelled 
in ID as decreasing a narrow potential barrier to zero height. In the adiabatic limit of 
this process, the phase step of = 2irrTvM/h 2 acquired by stirring reconnects at the 
nearest integer multiple of 2ir while the density is still low. In this way phase gradients 
are smoothed out before the density notch fills in completely leaving the condensate 
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in an azimuthally symmetric current-carrying state thus creating a pure toroidal ring 
current. In this idealized picture there is no upper limit for the current velocities. 
Above a critical velocity, however, the supercurrent is instable against phase slips on 
resting obstacles which might arise through imperfections in the trap |18|. 

It is not obvious that the simple ID picture of adiabatic stirring sketched above 
should be transferable to experimental conditions in 3D tori well out of the quasi- 
1D limit. The simulation of figure |l], however, shows that an essentially pure ring- 
current state can be generated by realistic choices of time scales and parameters. 
In order to get a rough estimate for the timescale of adiabaticity, consider the time 
£ a d = 27rrT/c max it takes for a sound wave to travel around the ring of radius rx once. 

We use the Bogoliubov speed of sound c max = \J Anh 2 aN pj M 2 psj] , which amounts 

to c max = 0.89mm s -1 at the peak density of the numerically obtained ground-state 
wavefunction, to estimate this time to i a d ~ 0.09s. It is further useful to compare the 
acceleration rate of the stirring paddle to a a d = c^ ax / '(27rr), which is the acceleration 
that brings the velocity up to the speed of sound in the time that it takes a sound wave 
to travel around the ring. We find that acceleration rates of about 0.1 a a d work well. 
Alternatively, instead of starting with the ground state of an intersected toroidal trap, 
the light sheet potential barrier can also be ramped up slowly. Our simulations show 
that ramping up the paddle linearly within about 3i a d is sufficiently slow in order to 
make no appreciable difference for the generation of ring currents. Further, retracting 
the core potential V^o, as realized in reference j|, leads to the generation of vortex 
filaments. We observe that essentially pure vortex states can be generated by ramping 
down the core-potential on a timescale of about 2i a d. 

A variation of the stirring experiment as described above can be used to generate 
solitons. In the non-adiabatic limit, when the intersecting light sheet is removed 
suddenly, the imposed phase step and density notch set an initial condition for 
the dynamics of the condensate. It has been discussed in other places, that the 
combination of a phase step [Q and a density notch is a very efficient way to start 
density- notch solitons p(| [2l|. For the example of a linear quasi- ID trap it is discussed 
by Carr et. al. [|TJ how the combination of a density notch generated by sudden 
removal of a Gaussian potential barrier in combination with a phase step can lead to 
the generation of one or more solitons with different velocities depending on the size of 
the phase step and the width of the potential barrier. We have performed simulations 
in tightly-confined toroidal traps with a transverse confinement in terms of healing 
lengths A/j between 2 and 5, where we observe a quasi-lD type of behaviour and the 
discussion in reference |2l|] applies. Our 3D simulations show the generation of band 
solitons upon sudden removal of the paddle. 

However, in situations where the confinement is less tight and the curvature of 
the torus becomes apparent the phenomenology changes: We observe the robust 
generation of a unique class of solitary excitations in a torus with transverse 
confinement in terms of healing lengths ~ 9. These excitations have properties 
of both solitons and vortices, as will be elaborated below, and therefore we call them 
solitonic vortices or svortices. Figure shows a simulated stirring experiment in a 
condensate of 10 6 atoms of Na where a phase step of approximately it is generated 
by stirring|J. Ramping down the paddle potential within 66ms (which is roughly 

|| We use the same potential as before with the following parameters for a condensate of 10 5 Na 
atoms: w x = 2tt21.0Hz, V c = 135nK, Ar c = 16.8/mi; V p = 50.6nK; Aa p = 3.75^m; Ar p = 30.0^m. 
The toroid radius in this simulation is = 36.0/mi. At peak density of the condensate ground 
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Figure 2. Stirring experiment where a svortex excitation is generated. Shown 
are theL density p and the phase of the condensate at different times as in 
figure hi. The first and the second column show the acceleration phase which 
lasts from t = to t = 0.66s. Afterwards the condensate is stirred with constant 
angular velocity ui = 27r0.032Hz. The paddle is contracted between t = 0.99s 
and t = 1.15s. Due to the transverse confinement being about 9 healing lengths, 
a dark soliton is not completely stabilized but transforms into a vortex by the 
influence of curvature. A state is generated which is on the transition between a 
dark soliton and vortex. 



1/3 tad with t a d ~ 0.23s in this simulation), a solitary wave is created that has 
characteristics of both solitons and vortices. While showing the phase singularity 
signature of a vortex, the density profile is squeezed by the tranverse confinement 
(see figure ||) . After the stirring potential is turned off, the svortex moves at constant 
angular velocity around the ring. Our simulations show that the velocity of the svortex 
is sensitive to the details of the stirring process like the paddle acceleration rate and 
timescale of paddle retraction. Similar to solitons, the velocity of svortices relative to 
the condensate is related to a phase offset. 

The velocity of solitons relative to the condensate is given by 

v = c max cos(A/2) (2) 

ji], [l6| H^], where A is the phase step across the density notch relating to standing 
black solitons (kink states) at A = it and to increasingly shallow gray solitons movin; 
with the speed of sound c max as A approaches zero. The experiments of reference 
have confirmed this relation for dark solitons in a BEC. We find that svortices obey a 
similar relationship. Figure^ shows the velocity-phase-offset relation for a svortex and 
a band soliton, which have been created by imaginary-time propagation in a long thin 
rectangular 2D box (of dimensions 8£ x 64£) with hard- wall boundary conditions . 
We find that the svortex velocity depends linearly on small deviations of the phase 
offset from tt resembling the relation (|J) for solitons. The proportionality coefficient, 
however, is below the value for solitons, given by the Bogoliubov sound speed c ma x- 
Simulations in differently sized boxes show that this coefficient depends strongly 
on the transverse confinement. We observe the discussed behaviour for transverse 
confinements Af^ between 6 and 12 healing lengths, ft is this velocity-phase-offset 

state in the torus we find a Bogoliubov sound speed of c max = 2.4mm s -1 and healing length of 
£ = 0.82/xm. 
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Figure 3. An unusual type of solitary excitations, which we call a svortex, 
can be generated in toroidal BECs by stirring. Shown are the iso-surface of the 
condensate at 16% of the maximum density and the phase <f> in the xy plane. 
The condensate wavefunction was prepared in a stirring experiment described in 
the text below. The solitary density deformation and phase singularity preserve 
their shape while rotating clockwise about the axis of the torus with roughly half 
the angular velocity of a ring current carrying a single quantum of vorticity. The 
current plot shows the wavefunction at t = 1.98s. 

relation, which makes the vortex confined to a thin waveguide a solitonic vortex 
or svortex. A detailed analysis of the svortex's properties including the collisional 
properties, which are considerably more complicated than those of solitons, will be 
published elsewhere p2| . We note that the term "vortex soliton" or "optical vortex 
soliton" has been used before in the nonlinear optics community in a much broader 
context pi whereas we introduce the term solitonic vortices referring to transversely 
confined vortices with specific solitonic properties. 

In conclusion, we have shown that ring currents and solitonic structures can be 
generated very efficiently in toroidal trap geometries by stirring the condensate with 
a blue detuned laser light sheet and consequently retracting the stirring paddle. Pure 
ring currents can be generated this way for a wide range of transverse confinement 
parameters. Density notch dark and gray solitons can be generated in quasi-lD 
transversely confined tori. In a regime of slightly wider confinement, a novel hybrid 
object with properties of solitons and vortices is generated. These solitonic vortices 
show a phase-offset-velocity relationship similar to solitons. 

We like to thank L. Carr and S. McKinney at the University of Washington and 
D. Feder and J. Denschlag at NIST for helpful discussions and comments. This work 
was partially supported by NSF Chemistry and Physics. J. B. acknowledges support 
from the Alexander von Humboldt Foundation through a Feodor Lynen Research 
Fellowship. 
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Figure 4. Velocity vs. phase-offset relation for a band soliton and a svortex in a 
2D box with dimensions 8£ X 64£. A stationary black band soliton (nodal-plane 
state) and a stationary svortex were set in motion by imprinting a phase step. 
The deviation of the measured phase offset A from 7r is approximately half the 
imprinted phase step. 
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